Electronic Structure and Correlation Effects in PuCoIns as Compared to PuCoGas 
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Since their discovery nearly a decade ago, plutonium-based superconductors have attracted considerable in- 
terest, which is now heightened by the latest discovery of superconductivity in PuCoIns . In the framework of 
density functional theory (DFT) within the generalized gradient approximation (GGA) together with dynami- 
cal mean-field theory (DMFT), we present a comparative study of the electronic structure of superconducting 
PuCoIns with an expanded unit cell volume relative to its PuCoGas cousin. Overall, a similar GGA-based 
electronic structure, including the density of states, energy dispersion, and Fermi surface topology, was found 
for both compounds. The GGA Pu 5/ band was narrower in PuCoIns than in PuCoGas due to the expanded 
lattice, resulting in an effective reduction of Kondo screening in the former system, as also shown by our DMFT 
calculations. 

PACS numbers: 74.25.Jb, 74.20.Pq, 71.27.+a, 71.28.+d 
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Introduction. The itinerant-to-localized crossover of the 
5/ electrons that occurs near plutonium in the actinide series 
is one of the most challenging issues in condensed matter. The 
cubic S phase of plutonium metal lies closer to the localized 
side of this boundary, similar to the heavy actinides (Am and 
beyond) where the 5/ electrons do not participate in bonding, 
while in a-Pu the 5/ electrons are itinerant and contribute to 
the bonding, similar to the light actinides (Th-Np) |[l|]. This 
change in bonding leads to a 25% larger volume in the S phase 
and a low-symmetry, monoclinic crystal structure for a-Pu, as 
well as a variety of unusual physical and mechanical proper- 
ties 10] ■ It is well established that these interesting phenomena 
in elemental Pu arise from the strong electronic correlation in 
the 5/ electrons fsl-Ml- With the discovery of superconduc- 
tivity in PuCoGas at Tc= 18.5 K and later in PuRhGas 
at Tc=^.l K inn , there is renewed interest in studying the 
strong electronic correlations that now also generate a transi- 
tion temperature an order of magnitude higher than in their 
CeMIns counterparts H, 12]. Furthermore, superconductivity 
(Tc ~ 2.5 K) has recently been discovered in PuCoIns 11311 . 
which has a unit cell volume 28% larger than its itinerant su- 
perconducting PuCoGas cousin, similar to the volume differ- 
ence between a-Pu and 5-Pu. Investigating these isostructural 
materials provides a particularly convenient way to probe the 
itinerant-to-localized crossover without the complication of a 
drastic structural change, and to help elucidate the origin of 
superconductivity in the Pu-based materials. 

In this Letter, we present a comparative study of the elec- 
tronic structure of PuCoIns and PuCoGas. Our calculations 
reveal they have the same number of Fermi surface sheets 
and similar band-center locations, although the details of their 
Fermi surface topology are sightly different. The expanded 
volume of PuCoIns causes a narrower bare 5/ band relative 
to that of PuCoGas. Moreover, the LDAh-DMFT calcula- 
tions show a reduction of Kondo screening in PuCoIns rel- 
ative to PuCoGas caused by the band narrowing, which tips 
the balance between competing Kondo and RKKY interac- 
tions towards magnetism and localization of the 5/ electrons 
in PuCoIns. 



Methodology. We performed electronic structure cal- 
culations of PuCoIns and PuCoGas within the framework 
of density functional theory (DFT) in the generalized gra- 
dient approximation (GGA) ll411 . Our calculations were 
carried out by using two relativistic band-structure meth- 
ods: The full-potential linearized augmented plane wave (FP- 
LAPW) method as implemented in the WIEN2k code 11511 . 
and the full-potential linear muffin-tin orbital (FP-LMTO) 
method as implemented in the RSPt code ITill . To address 
the 5/-electronic correlation issue, we used the GGA+U 
and GGAh-DMFT llTIl approximations, which are imple- 
mented in the WIEN2k code IH]. For the DMFT impurity 
solver, we used the vertex-corrected one-crossing approxima- 
tion (OCA) tH, which is reasonable for the description of 
more localized correlated electron systems. 

LDA bandstructure and Fermi surface topology. PuCoIns 
and PuCoGas crystalize in the tetragonal HoCoGas struc- 
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FIG. 1. (Color online) Calculated total energy versus volume for 
PuCoIns and PuCoGas in the paramagnetic state. The energy is 
shifted by 121004.4083 Rydberg for PuCoIns and 81612.8445 Ry- 
dberg for PuCoGas, respectively. The experimentally determined 
volumes are shown with arrows. 
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FIG. 2. (Color online) Calculated GGA total and partial density of 
states (DOS) for PuCoIns (top) and PuCoGas (bottom) in the para- 
magnetic state. The In 5p and Ga 4p DOS have been multiplied by a 
factor of 10 for clarity. 



ture (PA/mmm space group) with one internal z-coordinate 
for In or Ga. We calculated the volume dependence of the 
GGA-based total energy with the In and Ga z coordinates 
fixed at their experimental value of z(In)=0.306 ifisll and 
z(Ga)=0.312 iH, respectively. As shown in Fig. [l] we find 
the theoretical equilibrium volumes to be 1052.7 a.u.'^ for 
PuCoIns and 811.8 a.u.'^ for PuCoGas, which compare rea- 
sonably well with the experimental values of 1050.5 a.u."^ for 
PuCoIns and 820.2 a.u.'^ for PuCoGas. The good agreement 
at the GGA level indicates that the bonding between the tran- 
sition metal and ligand atoms is the dominant factor determin- 
ing the equilibrium volume of these Pu-1 15s, with the effect of 
the Pu 5/ electron correlation secondary in this regard, which 
is in striking contrast to the situation for elemental Pu jst). 
Hereafter, all calculations are performed at the experimentally 
determined lattice constants 11311 . 

Figure |2] shows the GGA total and partial density of states 
(DOS). Our results for the electronic structure of PuCoGas are 
in good agreement with earlier reports The two com- 

pounds exhibit somewhat similar features in the DOS. The 



strong spin-orbit coupling of Pu causes the 5/ states to be 
split into two manifolds or subshells, corresponding to a total 
angular momentum of j = 5/2 andj = 7/2. The partial DOS 
for Pu 5/ orbitals shows that the Pu 5/5 /2 states are the largest 
contribution at the Fermi energy, whereas the Co 3d and Ga 4p 
or In 5p orbitals have very small contributions. Furthermore, 
the narrow peak corresponding to Pu 5/5/2 is located slightly 
below the Fermi energy. Both the /2 and fj /2 peaks are nar- 
rower and exhibit less structure in PuCoIns than in PuCoGas, 
indicating a weakened hybridization in the former system due 
to the increased unit-cell volume. 

In Fig. [3] we show the band dispersion as a function of wave 
vector along high-symmetry lines. The Pu 5/ band character 
is indicated by the relative thickness of each line. The overall 
band structure of the two compounds is similar, and, as ex- 
pected from the DOS results, the bands in the vicinity of the 
Fermi energy consist mainly of Pu 5/ states. How these bands 
cut the Fermi energy determines the Fermi-surface topology. 
In total, there are four bands that cut the Fermi energy, which 
gives rise to four Fermi-surface sheets, as shown in Fig. |4] 
Among these four sheets, two of them are of hole character de- 
rived from the two lower bands cutting the Fermi energy and 
two of them of electron character derived from the two upper 
bands cutting the Fermi energy. Two hole pockets are centered 
at the r point in the zone center, while two electron pock- 
ets are centered at the AI point in the zone corners. At this 
level, the electronic structure of the Pu-1 15s bears some re- 
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FIG. 3. (Color online) Energy bands of PuCoIus (a) and PuCoGas 
(b). The thickness of the lines indicates the amount of Pu 5/ states 
present in each band. 




FIG. 4. (Color online) Calculated Fermi surface of PM PuCoIus (a) 
and FuCoGag (b). 
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TABLE 1. The Pu 5/ electron density within the muffin-tin sphere 
obtained in the GGA+f/ approximation for both SIC and AMF meth- 
ods for double counting corrections. 



semblance to the recently discovered Fe-based superconduc- 
tors ll23[l . Except for the small hole pocket, the Fermi surface 
exhibits a pronounced two-dimensional character, which is re- 
lated to the layered crystal structure with Pu forming a square 
lattice in each plane. A closer examination shows that the sec- 
ond hole Fermi surface (marked in red in Fig.|4|i, is less square 
on the z ~ ±7r zone face, which is due to the relative location 
of the four bands with respect to the Fermi energy (see Fig.|3]l. 
It is worth noting that the Fermi surfaces of the three known 
Pu-based superconductors are all qualitatively similar 

Pu 5f occupancy and correlations effects. Our GGA 
calculations for the Sommerfeld coefficient was found to be 
7GGA = 18 mJ/mol-K^ forPuCoIns and7GGA = 21 mJ/mol- 
for PuCoGas. These values are smaller by a factor of about 
10 and 5, respectively, than the experimental specific coeffi- 
cients, which are estimated to be 200 mJ/mol • ifisll for 
PuCoIngand 80 to liemJ/mol K^ |loil24l-|26ll for PuCoGa..^ . 
Although the renormaUzation effect is not as strong as in the 
Ce-115 compounds iIztI Esll . electronic correlations are still 
important. To understand how this affects the magnetism and 
superconductivity in the Pu-1 15s, it would be valuable to have 
some insight into the Pu 5/ valence of these compounds. For 
this purpose, we have performed GGA+U calculations by us- 
ing two different methods for double counting corrections; 
the self-interaction correction (SIC) approximation ll29ll and 
the around-mean field (AMF) method 13011 using an identical 
value of the muffin tin radius Rmt = 3.28 a.u. for both com- 
pounds. All calculations show that Pu 5/ weight remains at 



the Fermi level indicating some degree of mixed valent be- 
havior. The Fermi surface was qualitatively unchanged from 
that presented in Fig. 4 with the addition of U. Table I lists 
the 5/-orbitally projected electron density within the muffin 
tin sphere. Values close to 5 are consistent with previous es- 
timates based on the DFTh-DMFT method SH]. As can be 
seen, only relative occupations between compounds are mean- 
ingful, since the occupation depends on the basis sets used as 
well as on the double-counting coiTection method, and is sys- 
tematically larger with the AMF method than with the SIC ap- 
proximation. However, regardless of which scheme was used, 
was found to be larger in PuCoIns than in PuCoGas. A 



relatively larger value of ji^J indicates that the 5/ electron 
density is more localized in PuCoIns, because more 5/ elec- 
trons are pulled inside the muffin tin. 
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FIG. 5. (Color online) (a) Pu 5/ DOS at T = 20 K in the PM 
PuCoIns and PuCoGas from the GGA+DMFT calculations, (b) Ex- 
panded view near Ef- The data represented by blue line is for a 
hypothetical PuCoIns compound with a reduced unit cell volume by 
20%. 

The observation of an enhanced specific heat coefficient 
and a coherence feature in transport measurements indicate 
the importance of Kondo screening in these compounds, 
which is an effect that goes beyond what can be calculated 
in the framework of DFT within the local density based ap- 
proximation. In particular, the ultimate ground state of these 
compounds is determined by the competition between Kondo 
coupling and magnetic exchange interactions. To obtain a 
qualitative understanding of the Kondo exchange coupling in 
these systems, we performed GGAh-DMFT calculations. We 
used C/ = 4 eV for the Hartree component of the screened 
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Coulomb interaction, which is consistent with previous work 
on elemental Pu fsf-Sl- The remaining Slater integrals (F^, 
F^, and F^) were calculated using Cowan's atomic structure 
code fsil] and reduced by 30% to account for screening. Fig- 
ure |5] shows the Pu 5/ partial DOS, which exhibits a three- 
peak structure. The two broad peaks below and above the 
Fermi energy correspond to the j = 5/2 and j = 7/2 sub- 
shells, respectively, with an energy difference due mainly to 
the Hubbard U and the spin-orbit coupling. The central peak 
located very close to the Fermi energy is a Kondo resonance 
state, which is a hallmark of quantum many -body effects. This 
Kondo resonance, which constitutes a strongly renormalized 
quasiparticle band, is a generic feature that applies to both 
Pu-115 compounds. By comparing the renormalized band- 
width with the bare one, a rough estimate of the renormal- 
ization is about two orders of magnitude. This over-estimate 
is reasonable, since the impurity solver is based on the non- 
crossing type of approximation that underestimates the Kondo 
screening |17l] . Independent of the precise value, we clearly 
find a narrower quasiparticle bandwidth for PuCoIns than for 
PuCoGas. The fact that the quasiparticle band broadens to 
nearly the same amount as in PuCoGas, when the unit cell 
volume of PuCoIns is reduced by 20%, demonstrates that the 
reduction of Kondo screening is primarily caused by the ex- 
pansion of the lattice. 

The width of the renormalized band is a characteristic en- 
ergy scale that has been shown to control the maximum super- 
conducting transition temperature in the 115 materials [32 1 . 
Thus, the lower superconducting transition temperature of 
PuCoIns compared with PuCoGa5 may, in part, be a conse- 
quence of the reduction in Kondo screening. Future work will 
help elucidate the role of spin, orbital, and/or valence fluc- 
tuations for the observation of superconductivity in Pu-based 
compounds ll33ll34ll . 

Concluding remarks. We performed GGA band-structure 
calculations for the PuCoIns and PuCoGas superconductors. 
A similar electronic structure was found for both compounds. 
The expanded lattice in PuCoIns relative to PuCoGa5 results 
in a narrower bare Pu 5/ band width in PuCoIns and a conse- 
quent reduction in Kondo screenin g. W hen put in the context 
of the Doniach phase diagram i35ll36ll . our calculations sug- 
gest that PuCoIn5 is on the weak hybridization side, while 
PuCoGas is in the more strongly hybridized limit. We an- 
ticipate that a hypothetical PUC0TI5 compound would pos- 
sess a magneticly ordered state. To experimentally uncover 
the localization-delocalization transition of Pu 5/ electrons, 
PuCo(Ga,In)5 alloys would be natural candidates. This study 
supports the notion that an expansion in lattice constant can 
indeed drive the Pu 5/ electrons towards a localized state. 
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